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, dips in the transmission response. For all equal slits, the wavelengths of these sharp transmission 
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, minima can be varied by changing the width or the length of all the slits. In this paper we show 

Q ' that the width of the dip, i.e., the frequency range of minimum transmittance, can be controlled 

by making at least one slit different from the rest within a compound unit cell. In particular, we 

^ : 

investigate the effect that a change in the dielectric filling, or in the length of a single slit produces 



Phase resonances in transmission compound structures with subwavelength slits produce sharp 



in the transmission response. We also analyze the scan angle behavior of these structures by 
means of band diagrams, and compare them with previous results for all-equal slit structures. 
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I. INTRODUCTION 



The possibility of designing the transmission response of periodic structures with sub- 
wavelength slits has attracted the interest of the scientific community since these systems 
show promise for applications such as the characterization of attosecond pulses inte- 
grated polarizers |2| , optical data storage and external storage media. Also, a great effort 
has been devoted to the modeling of scattering by thick metallic plates with ID [s], [4, Is] 
and 2D 6l periodic perforations since the experimental demonstration of enhanced trans- 

PR 

mission by Ebbesen et al. [7|, 18|. Most of the studies carried out in connection with this 
phenomenon consider simple structures, i.e., periodic structures in which a unit cell is 
defined by a single hole or slit. 



Recent simulations reported in 



ll| suggest that compound gratings (periodic struc- 



tures comprising several slits in each unit-cell) can exhibit sharp transmission nulls when 
illuminated by p-polarized plane waves of particular frequencies. These sharp dips corre- 
spond to the excitation of phase resonances in the structure, which arise from a particular 
phase and amplitude distribution of the magnetic field in adjacent cavities of a unit cell, 
and are characterized by a field enhancement within the slits. The behaviour of phase res- 
onances in compound structures with subwavelength slits under oblique illumination has 
been studied recently HI], where the authors show that new phase modes, and therefore 
new transmission dips, that were not allowed for normal illumination due to symmetry 
reasons, appear for non-normal incidence. 



Phase resonances were 



irst reported in structures comprising a finite number of cavities 



on a perfect conductor 1^, Q, and later observed in compound reflection gratings under 



normal 



15l . Il6| and oblique [17|] illumination. Le Perchec et al. analyzed the excitation 
of these resonances in a two-slit system The additional parameters in the complex 
structure should allow for more degrees of freedom in the design of components such as 
nonlinear devices and frequency selective surfaces. Besides, for a selected region of the 
electromagnetic spectrum, the behavior of these complex structures can be compared to 
the effect produced by a defect in a photonic crystal (PC) [19]. Whereas a defect in 
a PC allows transmission within a forbidden band, a phase resonance in a compound 



transmission grating inhibits transmission in a highly transmitting region 



Recently, experimental demonstration of phase resonances in infinite compound gratings 



has been reported In th,s paper, the author, confirm earlier pred.et.ons made by 

Skigin and Depine in [9|] and show that a structure comprising three shts within each 
period, illuminated at normal incidence, exhibits a dip in the transmission response within 
each Fabry-Perot resonance peak. 

The theoretical approach used in this paper to model the ID compound transmission 
grating is based on the Generalized Scattering Matrix (GSM) formulation. It basically 
consists in expanding the fields in each region in its own eigenmodes and then matching 
them imposing the boundary conditions. Oblique excitations of both thick and thin sim- 
ple periodic frequency selective surfaces with arbitrary aperture cross-sections have been 



successfully analyzed with this approach 



21 



22 



23| . which has been recently extended 



to deal with periodic structures comprising complex unit cells [2J, |25(]. One advantage of 
the Extended Generalized Scattering Matrix (EGSM) formulation is that the individual 
GSM matrices describing any single discontinuity can be calculated separately and often 
a priori in a coordinate system best suited for an analytical solution. The assembly of 
GSMs describing a multitude of transverse discontinuities (i.e. multiple perforations per 
unit cell) forms a coordinate invariant EGSM that are easily cascaded via standard oper- 
ators for problems having arbitrary longitudinal complexity. This abstraction alleviates 
tedious coordinate notation and facilitates fast parallel computation of electromagnetic 
scattering from highly complex unit cells. 

The purpose of this paper is to show the design fiexibility that compound gratings give 
for getting particular responses. We investigate the effect produced by a change in the 
dielectric constant and in the length of a single slit within each period. We also analyze 
the dependence of the response on the angle of incidence. 

In Sec. II we summarize the extended generalized scattering matrix (EGSM) formulation 
used to solve the diffraction problem of a p-polarized plane wave by a compound struc- 
ture. In Sec. Ill we show numerical examples that illustrate the possibilities that these 
structures provide, by varying the dielectric filling as well as the length of the slits. We 
also show band diagrams that permit to visualize the angular dependence of the resonant 
frequencies and the bandwidth of the gaps, and compare them with those obtained for 
all-equal slits structures. Concluding remarks are given in Sec. IV. 



3 



II. THE EXTENDED GENERALIZED SCATTERING MATRIX FORMULA- 
TION 



The EGSM is an extension of the Mode Matching-Generahzed Scattering Matrix 
(MMGSM) method that allows computationally efficient solutions to a larger class of 
scattering problems. The method is capable of efficiently analyzing the reflection and 
transmission properties of TE- and TM- polarized plane waves incident at arbitrary angles 
onto thin or thick arbitrarily shaped frequency selective surfaces, gratings with generalized 
profiles, metal structures with multiple layers of dielectric fillings and/or coatings. The 
geometry of the scatterer can be periodic or aperiodic, transmitting and/or reflecting. 
In the present case, we consider an infinitely periodic one-dimensional compound struc- 
ture: each period comprises several perfectly conducting wires and three equal-width slits 
between wires are formed. Each slit can be filled with a different dielectric, as shown in 



Fig. 1. In what 
are found in [24 



bllows we summarize the fundamentals of the method. Further details 
25|. 

The structure is divided longitudinally into slices or waveguide sections. These sections 
and the interfaces between them are represented by their equivalent EGSM Sq and St, 
respectively, where the subscript denotes the number of a transition (T) or guiding section 
(G). To obtain the required transition and guided-section EGSMs, the fields in each slit 
are represented by modal expansions and appropriate boundary conditions are enforced 
at each transition: continuity of tangential electric fields in the open regions, null electric 
field on the perfect conductor, and continuity of tangential magnetic field. The equations 
resulting from these conditions are projected in convenient bases, which yields a matrix 
equation with the modal amplitudes towards and away from the transition as unknowns. 
Once all the guided sections and transition matrices are found, they are cascaded to find a 
relationship between the incident plane wave and the reflected and transmitted amplitudes 



24, 



25|. 



In the case of a uniformly filled waveguide with arbitrary longitudinal profile, where 
refiections from the gradual bends are negligible, it is sufficient to consider it as a straight 
section of the same waveguide having the length of the bent waveguide 26|, l27||. This is 
the case of the examples shown in Fig. 5. 
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III. RESULTS 



The results shown illustrate the evolution of the phase resonances in a perfectly conduct- 
ing compound thick-metal grating comprising three rectangular slits filled with different 
dielectrics. In the examples we consider symmetric structures in which the central slit 
can be different from the external ones in its dielectric filling (see Fig. 1) and eventually 
in its length (see inset of Fig. 5). A p-polarized plane wave impinges on the structure in 
all the examples shown. 

In Fig. 2 we show a sequence of curves when varying the dielectric constant of the central 
slit filling (Ec) while keeping all the rest of the parameters fixed: c/d = a/d = 0.08, 
h/d = 1.14, e = 1, 6q = 0°. We plot the zero-order transmittance as a function of the 
incident wavelength for six cases: Ec = 1, 1.2, 1.4, 1.6, 1.8, 2. The lower curve corresponds 
to identical slits, and the phase resonance appears as a sharp dip within the waveguide 
resonant peak of the slits. Phase resonances are produced in compound gratings due 
to the possibility of having different configurations of the field phase distribution within 



the different slits comprising each period of the structure [q]. In particular, for three 

slits within the period, only one resonance is expected for normal incidence: the (H h) 

mode, i.e., the mode in which the phases of adjacent slits are reversed. The number of 
possible phase configurations is finite and depends on the number of slits. A mathematical 



criterion governing the resonances using the EGSM formalism is given in 2m, l25 |. 
In general, two requirements are needed for a phase resonance to occur: i) at least one 
non-zero phase difference is found between the field phases in adjacent slits — what is not 
allowed in simple periodic structures — and ii) a particular distribution of the field ampli- 
tude, naturally generated by the incidence conditions, is obtained [ij] . Phase resonances 
are also characterized by a strong intensification of the field inside the slits or inside the 



cavities, in the case of compound reflection gratings 



3- 



In particular, when the phases 



in adjacent slits are opposite to each other, vr resonances can be excited Regarding 
each slit as a waveguide, these two conditions are usually fulfilled within the waveguide 
mode resonant wavelength. 

It can be observed that as the permittivity of the central slit is increased, the width 
of each phase resonance dip also increases. This suggests that the bandwidth of the 
narrow gaps generated by means of phase resonance excitations can be controlled. This 



characteristic makes this structure attractive from the point of view of several apphcations, 
such as filters and polarizers. It can also be regarded as a complementary structure to 
the photonic crystal. A typical band diagram of a photonic crystal presents band gaps, 
i.e., frequency ranges that are not allowed to propagate inside the crystal and therefore, if 
the structure is illuminated by a plane wave of a frequency within the gap, all the power 
is reflected 191]. However, if the perfect periodicity of the photonic crystal is broken by a 
defect in the structure, allowed states arise within the gap, enabling transmission within 
the originally forbidden gap. In the compound gratings considered in this paper, the 
structures are essentially transmitting, at least within the waveguide resonance peaks, 
as observed in Figure 2. Thus, when a simple grating is illuminated by a plane wave 
of a wavelength within the waveguide resonance peak, most of the power is transmitted 
through the structure. The addition of complexity to the period in the form of a group 
of slits can be regarded as a defect in the perfect periodicity, and in this case one or more 
forbidden channels are formed within the allowed bands 111 ]. 

As it was already investigated in for all-alike slits structures, when compound gratings 
are obliquely illuminated, new possiblities of phase configurations within the slits open 
up, resulting in new phase resonances. For instance, for three equal slits we have the 

(H ) and the (+ + -) modes in addition to the one we had for normal incidence. To 

visualize the effect of changing the dielectric filling of the central slit in the response of 
the structure under oblique illumination, we show in Fig. 3 the transmittance for the 
same structure considered in Fig. 2, but for an oblique incidence: 6*0 = 30°. It can 
be observed that now there are two dips within each waveguide resonance maximum, 
i.e., around A/rf ~ 1.5 and 2.5. For oblique incidence, the symmetry condition imposed 
by normal illumination is removed, and this allows new phase configurations inside the 
slits, which produce new dips in the transmission response. The resonance dip at \/d 
slightly smaller than 2.5 for all equally filled slits (lower curve in Fig. 3(a)) shifts to 
larger wavelengths with decreasing quality factor [Q) for small variations of the central 
dielectric constant Sc (see the curve for = 1.1). From then on, it collapses into a very 
high-Q resonance for X/d slightly larger than 2.5, and stays fixed throughout all further 
increases of Ec- The dip for X/d slightly larger than 2.5 in the lower curve {Ec = e) shifts 
to larger wavelengths with decreasing Q as Ec increases. This behaviour suggests that 
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both dips behave differently for variations in Ec- The persistence of the high-Q resonance 
dip (leftmost dip within the peak) for all dielectric variations suggests that it is fixed 
by geometry and excitation, whereas the Q and the position of the rightmost dip are 
dielectric filling dependent. Besides, the transmittance peak to the left of the sharp null 
does not change much as £c is varied, whereas the peak to the right of the rightmost null 
is certainly tunable by changing the relative dielectric constant of the central slit filling. 
To get more insight about the influence of changing the central slit filling when the 
structure is obliquely illuminated, in Fig. 4 we show contour plots of the zero-order trans- 
mittance as a function of ctod and ujd/c, where ckq is the incident wave vector component 
along the direction of periodicity and uj is the frequency of the incident wave. The lighter 
zones represent higher transmitted intensities. Fig. 4(a) corresponds to a compound 
structure with all-alike slits. It can be observed that there are mainly two transmission 
peaks, for uod/c ~ 2.5 and 5. These maxima are scored by nulls, the number of which 
depends on the incidence angles: for normal illumination there is one dip within each 
peak, whereas for larger angles two dips are found within each peak. When the dielectric 
constant of the central slit filling is increased in steps of 0.2, we get the plots of Figs. 
4(b) to 4(f) [Sc = 1.2, 1.4, 1.6, 1.8 and 2). It can be noticed that for normal incidence, 
the dip width increases with £c, as it was already observed in Fig. 2, and this trend is 
also maintained for the larger wavelenght dip at oblique incidence, as observed in Fig. 3 
for Oq — 30°. The results shown in Fig. 4 suggest that the phase resonance mechanism 
could be exploited not only to design structures with a prescribed number of resonances 
at specified wavelengths, but also to control their bandwidth by changing the dielectric 
filling. 

The effect produced by changing the dielectric filling of the central sht in the transmitted 
response of the structure is similar to what occurs when the length of the central slit is 
changed. The propagation constant of each mode of a slit bounded by perfectly conducting 
walls and filled with a material with dielectric constant e\ is: 

where C is the speed of light in vacuum and m is an integer. Then, for a given length of the 
slit Izi, the phase difference gained by the mode when going through the slit is Izi. If 
the slit is now filled with another dielectric with Si — olSx^ the new propagation constant 
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is 

and then the same phase difference would be obtained for another slit length lz2 such that 

j3i^ IZi = /32m lZ2 • (3) 

In particular, if the slit has a subwavelength width, as in the present case, the only 
propagating mode is the first mode, which corresponds to m = in ([1]) and ([2]). Then, 
condition ([3]) yields Izi = y/a\z2. Taking into account the above relationship between the 
length and the dielectric constant of the slit filling, it is to be expected that the response 
of a structure with a higher central permittivity should be similar to that of the structure 
comprising equally filled slits, but with a longer central slit. The EGSM method used to 
solve the diffraction problem allows for the treatment of bent waveguides with arbitrary 
longitudinal profile. Then, we increased the length of the central slit (Iz) by bending it 
slightly, as shown in the inset of Fig. 5. The results in Fig. 5 show the evolution of the 
transmitted intensity as the central slit length is increased, when keeping e = Ec = 1, 
for the same parameters considered in Fig. 2. It can be observed that the three curves 
(Iz = h/d, \/1.2h/d and y/lAh/d) for Ec = 1 are very similar to those of Fig. 2(a), for 
Ec = 1, 1.2 and 1.4 and Iz = h/d. This result also shows that the resonance bandwidth 
can be controlled not only by changing the dielectric filling, but also by varying the length 
of the slits. 



IV. CONCLUSIONS 



We have investigated the response of compound transmission wire gratings with subwave- 
length slits under oblique incidence. To solve the diffraction problem we used the EGSM 
method, which is very versatile. The ease of its formulation significantly reduces the time 
and effort required for more advanced structures. We focused in a structure comprising 
three slits in each period and analyzed the influence that a change in the dielectric filling 
and in the length of the central slit has over the resonant dip. We have shown numerical 
results that evidence the possibility of controlling the bandwidth of phase resonance dips 
by properly choosing the geometrical and constitutive parameters. The angular behavior 
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of these structures was also analyzed by means of band diagrams, and the results obtained 
were compared with previous ones corresponding to all-equal slit structures. The capabil- 
ity of exciting and controlling the bandwidth of phase resonances in compound gratings 
opens up new possibilities for practical applications, such as polarization sensitive aper- 
ture shapes for field enhancement devices. 
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FIG. 1: Sketch of a three-slit compound grating. The relative permittivity of the material inside 
the middle slit differs from that in the edge slits. In this work, the slit length, h, and the slit 
widths, c and a, are varied and the resulting response is found as a function of the incidence 
angle ^o- 
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Xld 

(b) 

FIG. 2: Zero-order transmittance as a function of the normalized incident wavelength for a 
structure with a three slit unit cell, incidence angle do = 0°, for six values of the dielectric 
constant of the central slit filling ec=l, 1-2, 1.4, 1.6, 1.8, 2. The rest of the parameters are: 
c/d = a/d = 0.08, h/d = 1.14 and e = 1. 
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(b) 

FIG. 3: Zero-order transmittance as a function of the normalized incident wavelength for a 
structure with a three slit unit cell, incidence angle = 30°, for several values of the dielectric 
constant of the central slit filling £c='\-, 1-1, 1-15, 1.2, 1.4, 1.6, 1.8, 2. The rest of the parameters 
are: c/d = a/d = 0.08, h/d = 1.14 and e = 1, where d is the period. 



14 



(a) 



(b) 





FIG. 4: Contour plots of the zero-order transmitted intensity as a function of qq d, and to d/c, 
for the same cases considered in Fig. 2: (a) Sc = 1; (b) Ec = 1.2; (c) £c = 1.4; (d) Ec = 1.6; 
(e) £c = 1-8; (f) £c = 2. Black represents null transmittance and white represents total (unity) 
transmittance. 
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FIG. 5: Zero-order transmittance as a function of the incident wavelength for a structure with 
a three slit unit cell, incidence angle 6q = 0°, for three values of the central slit length Iz = h/d, 
V^h/d, ^/lAh/d. The rest of the parameters are: c/d = a/d = 0.08, h/d = 1.14 and £ = £c = 
1. The inset shows the scheme of a structure with a longer central slit. 
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